Limits on spin-independent WIMP-nucleon interactions from the two-tower run of the 

Cryogenic Dark Matter Search 
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We report new results from the Cryogenic Dark Matter Search (CDMS II) at the Soudan Under- 
ground Laboratory. Two towers, each consisting of six detectors, were operated for 74.5 live days, 
giving spectrum- weighted exposures of 34kg-d for germanium and 12kg-d for silicon targets after 
cuts, averaged over recoil energies 10-100 keV for a WIMP mass of 60 GeV/c 2 . A blind analysis was 
conducted, incorporating improved techniques for rejecting surface events. No WIMP signal exceed- 
ing expected backgrounds was observed. When combined with our previous results from Soudan, the 
90% C.L. upper limit on the spin-independent WIMP-nucleon cross section is 1.6xl0 -43 cm 2 from 
Ge, and 3xl0~ 42 cm 2 from Si, for a WIMP mass of 60GeV/c 2 . The combined limit from Ge (Si) 
is a factor of 2.5 (10) lower than our previous results, and constrains predictions of supersymmetric 
models. 

PACS numbers: 14.80.Ly, 95.35.+d, 95.30.Cq, 95.30.-k, 85.25.Oj, 29.40.Wk 



One quarter of the energy density of the universe con- 
sists of non-baryonic dark matter 0, which is gravita- 
tionally clustered in halos surrounding visible galaxies Q • 
The Weakly Interacting Massive Particle (WIMP) Q, a 
dark matter candidate, arises independently from con- 
siderations of Big Bang cosmology and from supersym- 
metric phenomenology, where the neutralino can be a 
WIMP • A WIMP has a scattering cross section with 
an atomic nucleus characteristic of the weak interaction 
and a mass comparable to that of an atomic nucleus. 

The CDMS II experiment 0, Q is designed to detect 
recoils of atomic nuclei that have been scattered by in- 
cident WIMPs in germanium (Ge) and silicon (Si) crys- 
tals. Events with recoil energies of a few tens of keV 
and rates <1 event kg -1 d _1 are expected 0, The 
CDMS II search is most sensitive to spin-independent 
(SI) WIMP-nucleon scattering amplitudes Coherent 
superposition of SI amplitudes gives Ge better sensitivity 
than Si, except for small WIMP masses, where kinemat- 
ics increase the Si sensitivity. In this report we interpret 
our data with the SI ansatz. Another report describes a 



spin-dependent interpretation of our data |lfj. 

The CDMS II apparatus in the Soudan Underground 
Laboratory has been described elsewhere [|| U At 
the experiment's core, Z(depth)-sensitive Ionization and 
Phonon detectors (ZIPs) measure the ionization and 
athermal phonon signals caused by recoiling particles in 
Ge and Si crystals [3, EI • We report new results from 
the most recent CDMS data run collected between March 
25 and August 8, 2004, consisting of 74.5 live days. Six 
Ge (250 g each) and six Si (100 g each) ZIPs were oper- 
ated in two vertical stacks ( "towers" ) at a temperature of 
50 mK. This report includes the first data from Tower 2, 
which contains four Si and two Ge ZIPs. Improvements 
made since the previous run [|| include a reduction of 
electrical noise and more frequent infrared illumination 
to clear the crystals of space-trapped charge. 

The ZIP detector provides event-by-event discrimina- 
tion of nuclear recoils from the dominant background of 
electron recoils. The ratio of ionization energy to phonon 
energy ("ionization yield") is ~ 0.3 for nuclear recoils, 
normalized to electron recoils (see yields 1 in Fig. 1). 
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Electron recoils near the detector surface suffer from poor 
ionization collection and can mimic the ionization yield 
of nuclear recoils that occur throughout the detector. 
Each ZIP's phonon sensors are divided into four quad- 
rants. Timing and signal amplitude comparisons among 
the quadrants provide discrimination against electron re- 
coils, particularly those near the surface. Also, surface 
electron recoils often deposit energy in adjacent ZIPs 
within a tower, causing multiple-detector events. Energy 
from a WIMP recoil would be contained in one detector. 

To avoid bias, we performed blind analyses. Events 
in and near the signal region were removed from WIMP- 
search data sets, or "masked" 01- The cuts that define a 
signal were determined using calibration data from 133 Ba 
and 252 Cf sources and from the non-masked WIMP- 
search data. Seven million electron recoils were collected 
using a 133 Ba source of gamma rays, exceeding the com- 
parable WIMP-search data by a factor of twenty. Half of 
the 133 Ba data were used to define analysis cuts and the 
other half to test cut definitions and estimate expected 
backgrounds. The detector response to nuclear recoils 
was characterized with 10 4 neutrons from a 252 Cf source, 
collected in four separate, several-hour periods during the 
run. 

Data from the 133 Ba source were used to monitor de- 
tector stability and to characterize detector performance. 
One Ge detector in Tower 2, ZIP 5 (T2Z5(Ge)), had a 
spatial region of abnormal ionization response which we 
excluded from analysis. The Si detector T1Z6, known to 
be contaminated with 14 C, a beta emitter, was excluded, 
as were detectors TlZl(Ge) and T2Zl(Si) due to poor 
phonon sensor performance. We report results from the 
remaining 5 Ge and 4 Si ZIPs, chosen before unmasking 
the WIMP signal region. 

Most event selection criteria, based on quantities from 
the four phonon and two ionization pulses from each de- 
tector,are similar to those described in our previous re- 
ports 0,0} ■ However, we completed five distinct analyses 
|l3j focused on improving existing methods and devel- 
oping new techniques to reject surface electron recoils. 
Events with low ionization yield in the 133 Ba calibration 
data, which are from surface electron recoils, were used 
to develop rejection criteria. All cuts were frozen prior 
to unmasking the signal region. 

Phonon pulses from surface recoils are more prompt 
than those from recoils in the detector bulk. Two tim- 
ing quantities in the quadrant with the largest phonon 
signal or "local quadrant" are particularly powerful: the 
time delay of the phonon signal relative to the fast ion- 
ization signal, and the phonon pulse risetime 0, Q] . For 
the first and simplest of the five analysis techniques, the 
sum of delay and risetime forms a timing parameter, af- 
ter energy corrections to delay and risetime are applied 
to achieve an energy-independent distribution. Figure 
shows the ionization yield versus this timing parameter 
for a typical Ge detector. The 133 Ba calibration events 
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FIG. 1: Ionization yield versus timing parameter (see text) 
for calibration data in T2Z3(Ge), with recoil energies in 
the range 10-100 keV. Typical bulk-electron recoils from the 
133 Ba source of gamma rays are dots (red), with yield near 
unity. Low-yield 133 Ba events (+, black), attributed to sur- 
face electron recoils, have small values of timing parameter, 
allowing discrimination from neutron-induced nuclear recoils 
from 252 Cf(o, blue), which show a wide range of timing pa- 
rameter values. The vertical dashed line shows the minimum 
timing parameter allowed for WIMP candidates, and the box 
shows the approximate signal region, which is in fact weakly 
energy dependent. (Color online.) 

with yields of 0.1 to 0.8, identified as surface electron re- 
coils, show a smaller timing parameter than most nuclear 
recoils induced by the 252 Cf source. 

For data from the Ge detectors, we require that can- 
didates for WIMP-induced nuclear recoil exceed a mini- 
mum value for this timing parameter (see Fig. [2] upper- 
right). Because this "timing parameter" analysis is sim- 
ple and robust, we agreed before unmasking to report its 
results. The expected sensitivities computed for the four 
more advanced analyses described below were all within 
±20% of that computed for this technique. 

Two of the four advanced analysis techniques evaluate 
the goodness of fit, x 2 , for surface electron versus bulk 
nuclear recoil hypotheses. These methods use three vari- 
ables including their correlations: time delay, risetime, 
and "partition" . The partition parameter, a measure of 
energy distribution between the phonon quadrants, is de- 
fined as the ratio of phonon energy in the local quadrant 
to that in the diagonal quadrant. The difference in 
between surface electron and bulk nuclear recoil hypothe- 
ses, A^ 2 , is used as the principal discrimination parame- 
ter in these analyses. One method uses energy-dependent 
variances and the other energy-independent variances. 

For the Si detectors, we decided after unmasking to use 
the energy-dependent \ 2 technique, which has the best 
expected sensitivity to a nuclear recoil signal from low- 
mass WIMPs for any of our five methods. This analysis 
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technique is four times as sensitive as the timing param- 
eter analysis for low-mass WIMPs. Before unmasking, 
we set the minimum requirement on the A\ 2 for surface 
electron recoil rejection (see Fig. flower-left). These cri- 
teria are used for our WIMP-search results for Si. 

The two remaining analysis techniques confirm the ro- 
bustness of our Ge and Si results. One technique com- 
bines delay, risetime and partition in a neural net analy- 
sis, and the other technique exploits additional informa- 
tion from the fitted signal pulses to reconstruct recoil po- 
sition and type. These two and the \ 2 techniques promise 
further improvements in sensitivity for the larger expo- 
sures planned in our future runs, beyond the improve- 
ments already demonstrated in Fig. 
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FIG. 2: Variables used to reject surface electron recoils (7- 
100 keV), for data from T2Z4 (Si/left) and T1Z2 (Ge/right), 
for the timing parameter (upper) and A\ 2 (lower). Light 
(red) lines show distributions of low-yield electron recoils from 
the 133 Ba source, while dark (blue) lines show distributions of 
nuclear recoils from the 252 Cf source. Dashed lines indicate 
the minimum values for acceptable WIMP candidates. A cut 
on the timing parameter is used for the Ge detectors, while 
a requirement on A^ 2 is used for the Si detectors. (Color 
online.) 



Upon unmasking the Ge data, one candidate with 
10.5 keV recoil energy was found to pass all cuts in the 
timing parameter analysis. Unmasking the Si data re- 
vealed that no events passed all cuts. Figure shows the 
unmasked data in ionization yield versus recoil energy; 
the data before and after application of the cut intended 
to reject surface electron recoils are shown. All analysis 
techniques showed consistent results. 

After unmasking the Ge data, we realized that the can- 
didate occurred in a detector during an interval of time 
when that detector suffered inefficient ionization collec- 
tion. This defect by itself would prevent us from claiming 
the event was evidence for a WIMP-induced nuclear re- 
coil. The candidate is also consistent with the rate of 
expected background. Although we do not claim this 
event as a nuclear recoil from a WIMP, we do include it 
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FIG. 3: Ionization yield versus recoil energy for events in all 
Ge detectors (upper) and all Si detectors (lower) passing ini- 
tial data selection cuts prior to applying the surface electron 
recoil rejection cut. The signal region consists of recoil ener- 
gies exceeding 7keV, shown with a vertical dashed line, and 
yields between the curved lines defined with recoils induced by 
the 252 Cf source while WIMP-search data were still masked. 
Below 7keV, separation between nuclear and electron recoils 
becomes poor. Events passing the surface electron recoil cut 
are a star (red) inside the signal region and dark filled circles 
(blue) outside. Bulk-electron recoils with yield near unity are 
above the vertical scale limits. (Color online.) 

in setting limits on the WIMP-nucleon cross section. 

The expected surface-event backgrounds are estimated 
by multiplying two factors. The first factor is the num- 
ber of events in the signal region before surface-event 
cuts as obtained upon unmasking. The second factor 
is the fraction of surface events expected to pass the 
surface-event cut. This fraction may be estimated with 
the actual passing fraction of low-yield events similar to 
the single-detector event background. An estimate using 
the passing fraction from 133 Ba Ge (Si) data indicates a 
surface-event background of 0.13±0.05 (0.90±0.4) events. 
However, we decided before unmasking to base our 
background estimate on the passing fraction of WIMP- 
search multiple-detector events in a wide low-yield re- 
gion. The number of surface events expected to pass the 
surface-event cuts is 0.4±0.2 (stat) ±0.2 (sys) between 
10-100 keV in Ge and 1.2±0.6 (stat) ±0.2 (sys) between 
7-100 keV in Si. Evidence suggests that beta decays 
of radioactive nuclides distributed across the detectors 
cause most surface electron recoils in the WIMP-search 
data. The spatial distribution of these contaminants dif- 
fers from that of surface electron recoils from the 133 Ba 
source, which might explain the difference in estimated 
background levels. From simulation methods reported in 
, the expected background due to cosmogenic neutrons 
that escape our muon veto is 0.06 events in Ge and 0.05 
events in Si. 
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merit of Energy. The ZIP detectors were fabricated in the 
Stanford Nanofabrication Facility operated under NSF. 
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FIG. 4: WIMP-nucleon cross section upper limits (90% C.L.) 
versus WIMP mass. The upper CDMS Ge curve also uses 
data from the current run, while the lower Ge curve in- 
cludes the previous run 6). Supersymmetric models allow 
the largest shaded (light-blue) region [l5|. and the smaller 
shaded (green) region [TH ]. The shaded region in the upper 
left (see text) is from DAMA 0, and experimental limits 
are from DAMA O, EDELWEISS [H, and ZEPLIN (U- 
(Color online.) 

Figure 0] shows the upper limits on WIMP-nucleon 
cross sections calculated from the Ge and Si analyses 
reported here using standard assumptions for the galac- 
tic halo ||. For the upper Ge limit, data between 10- 
100 keV from this run are used. Also shown is the com- 
bined limit obtained from this report and our earlier work 
0,0- For the combined Ge limit, we have included data 
in the 7-10 keV interval of recoil energy from the run 
reported here |21| . The combined result for Ge limits 
the WIMP-nucleon cross-section to < 1.6 x 10~ 43 cm 2 
at the 90% C.L. at a WIMP mass of 60GeV/c 2 , a fac- 
tor of 2.5 below our previously published limits. This 
new Ge limit constrains some minimal supersymmet- 
ric (MSSM) parameter space and for the first time 
excludes some parameter space relevant to constrained 
models (CMSSM) 0. 

The Si limit in Fig. 0] is based on standard halo as- 
sumptions using Si data from 7-100 keV in this run. 
The Si result limits the WIMP-nucleon cross-section to 
<3xl0~ 42 cm 2 at the 90% C.L. at a WIMP mass of 
60GeV/c 2 . This Si result excludes new parameter space 
for low-mass WIMPs, including a region compatible with 
interpretation of the DAMA signal (2-6 and 6-14keVee 
bins) as scattering on Na (lif. 
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